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Abstract 

Ultra high energy neutrinos are produced by the interaction of hadronic cosmic 
rays with the cosmic radiation background. More exotic scenarios like topological 
defects or new hadrons predict even larger fluxes. In particular, Earth-skimming 
tau neutrinos could be detected by the Fluorescence Detector (FD) of Pierre Auger 
Observatory. A detailed evaluation of the expected number of events has been 
performed for a wide class of neutrino flux models. An updated computation of the 
neutrino-nucleon cross section and of the tau energy losses has been carried out. 
For the most optimistic theoretical models, about one Earth-skimming neutrino 
event is expected in several years at FD. 



PACS numbers: 95.85.Ry, 13.15. +g, 96.40. Tv, 95.55.Vj, 13.35.Dx; 



1 Introduction 



The measurement of Ultra High Energy Cosmic Rays (UHECR) flux is the goal of a 
wide class of past, present and future detectors • UHE neutrinos are expected to be 
produced by the interaction of hadronic matter with the surrounding radiation/matter. 
A search for this signal is currently performed by several Neutrino Telescopes [TU]-|16j. 

Neutrinos with energy above 10 17 eV are expected to originate from the interac- 
tion of UHE cosmic rays with the Cosmic Microwave Background (CMB) via the 7r- 
photoproduction, p + jcmb — ► n + n + , the so-called cosmogenic neutrinos [Ej. The 
prediction for such a flux is however affected by several uncertain physical quantities, 
namely the spatial distribution of astrophysical sources, the ejected proton fluxes (if pro- 
ton) and the way of modelling the diffuse extragalactic electromagnetic background in 
the different frequency regions. One can assume a reasonable ansatz for all these quan- 
tities combined with the measurement of the diffuse photon flux in the GeV region by 
EGRET [18,, and the AGASA/HiRes data. These models and the predictions of more 
exotic scenarios have been exhaustively discussed in several papers (see for example Ref.s 

mm)- 

High energy neutrinos are hardly detected, as they are almost completely shadowed 
by Earth and rarely interact with the atmosphere. An EeV neutrino has an interaction 
length of the order of 500 km water equivalent in rock and, even crossing horizontally the 
atmosphere (360 meters water equivalent), only one neutrino out of thousand will be in- 
teracting. Due to the very low expected flux and the small neutrino-nucleon cross section, 
km 3 -neutrino telescopes and giant surface arrays have very few chances of detection. 

In this framework, an interesting strategy for v T detection is described in Ref.s |21j - 
[3*2] . As shown for example in Figure 1 of Ref. [2H] 5 for energy between 10 18 and 10 21 eV 
the r decay length is not much larger than the corresponding interaction range. Thus, 
an energetic r, produced by Charged Current (CC) v T interaction not too deep under the 
surface of the Earth, has a chance to emerge in the atmosphere as an upgoing particle. 
Unlike r's, muons crossing the rock rapidly loose energy and decay. Almost horizontal 
v T , just skimming the Earth surface, will cross an amount of rock of the order of their 
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interaction length and thus will be able to produce a corresponding r, which might shower 
in the atmosphere and be detected. In order to estimate the number of upgoing r expected 
in the Pierre Auger Observatory (PAO), one needs to know the value of the neutrino- 
nucleon cross section for CC interaction, u T + N — > r + X. 

The aim of this paper is to estimate the number of possible upgoing r showers which 
the Fluorescence Detector (FD) of PAO could detect. The predictions are analyzed with 
respect to their dependence on different neutrino fluxes and by using a new estimate of 
v — N cross section. 

The paper is organized as follows. In section 121 different models for neutrino flux 
predictions are discussed. In section El the general features of deep inelastic neutrino 
cross sections are outlined for neutrino energies up to 10 21 eV. In section HJ high energy 
r propagation through matter is illustrated by considering all the relevant interaction 
mechanisms. Average values for r energy loss are provided by taking into account the 
recent calculations of photonuclear interaction given in Ref. [33]. In section the number 
of expected upgoing r showers is derived and discussed. Finally in section El we give our 
conclusions and remarks. 

2 Neutrino flux estimates 

Ultra High Energy protons, with energy above ~ 10 20 eV, travelling through the universe 
mostly loose their energy via the interaction with CMB radiation. The large amount of 
charged and neutral pions produced will eventually decay in charged leptons, neutrinos 
(cosmogenic neutrinos JEj), and high energy gamma rays. 

At the GeV energy range the extragalactic diffuse gamma-ray background was mea- 
sured by the EGRET experiment |18j . This measurement provides an upper bound for 
possible neutrino fluxes from pion production. In particular, it gives the expected maxi- 
mum flux of cosmogenic neutrinos from an initial spectrum of measured UHE protons [35 . 
It is worth noticing that, since at least part of UHECR are protons, the existence of cos- 
mogenic neutrinos is guaranteed, even if their flux is very uncertain. In Figure Qthe GZK 
neutrino flux for three possible scenarios is plotted. The thick solid line gives the case of 
an initial proton flux oc 1/E, by assuming in addition that the EGRET flux is entirely 
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Figure 1: Cosmogenic neutrino fluxes as a function of energy. Thick solid line is for 
an initial proton flux oc l/E, by assuming that the EGRET flux is entirely due to tt- 
photoproduction (GZK-H). Thin solid line shows the neutrino flux when the associated 
photons contribute only up to 20% in the EGRET flux (GZK-L). Dashed line stands for 
an initial proton flux oc l/E 2 (GZK-WB). The experimental points represent the UHECR 
flux measured by AGASA (crosses) and HiRes (stars). 



due to 7r-photoproduction (GZK-H). The thin solid line shows the neutrino flux when the 
associated photons contribute only up to 20% in the EGRET flux (GZK-L). The dashed 
line stands for the conservative scenario of an initial proton flux oc l/E 2 (GZK-WB). In 
this case the neutrino flux is compatible with the so-called Waxman-Bahcall limit [36,. 
Note that no lower bounds can be set for the cosmogenic neutrino flux. In particular, 
in the most conservative but rather unrealistic case, the astrophysical sources cannot ac- 
celerate protons up to energies above GZK cutoff, and thus the secondary neutrinos will 
be produced in negligible quantities. All neutrino fluxes presented in this section were 
calculated by a propagation code [37] which takes into account the neutrino production 
via 7r-photoproduction with microwave, infrared, optical and radio photon backgrounds, 
as well as in neutron decay. 

Most of the models trying to explain highest energy cosmic rays (E > 10 20 eV) in terms 
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Figure 2: Neutrino fluxes in exotic UHECR models. Solid line represents the neutrino 
flux prediction in a model with new hadrons (NH) [SB]; whereas the dashed line is the 
neutrino flux for a topological defect model (TD) [39 . 



of exotic particles, predict a large associated flux of neutrinos. In Figure El the expected 
neutrino flux for two of such scenarios is plotted. One of them is the model of new hadrons 
(NH) [3H] , with mass M ~ 2 — 5 GeV, capable of generating UHECR events above GZK 
cutoff. In SUSY theories, for example, the new hadrons are bound states of light bottom 
squarks or gluinos and, once produced in suitable astrophysical environments, can reach 
the Earth without significant energy losses. In spite the production of new hadrons is a 
subdominant process, it generates a large number of neutrinos (see Figure |2J). 

The dashed line in Figure El shows the neutrino flux for a Topological Defects model 
(TD) (for a review see [39J). In this case UHECR events with energy E > 10 20 eV are 
explained in terms of 7's which are produced in the decay of heavy particles with mass 
of the order of 10 22-23 eV. As in the previous case, the associated neutrino flux for this 
kind of models is extremely large. 

We do not discuss here the so-called Z-burst scenarios, which attempts to explain 
UHECR as products of Z-boson decay. These models are in fact strongly disfavored [20] 
by the upper bounds on UHE neutrino flux put by FORTE [lUj an d GLUE jJT] and by 
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the cosmological limits on neutrino mass set by WMAP and LSS data [J21II31E]- 

In the following sections we will estimate the sensitivity of PAO to the UHE tau 
neutrino flux both in the case of cosmogenic neutrinos (Figure [Q) and of exotic models 
(Figure EI). 

3 Neutrino-Nucleon cross section in the extremely 
high energy limit 

At energy above 1 GeV neutrino-atoms interaction is dominated by the process of Deep 
Inelastic Scattering (DIS) on nucleons, since the contributions of both elastic and quasi- 
elastic interactions become negligible. The effect of the neutrino scattering with atomic 
electrons will not be taken into account here, since the cross section for this process is, 
at each energy, about three orders of magnitude lower than the neutrino-nucleon cross 
section 1 . 

Detectable leptons are produced through Charged Current interaction, 

^) + JV-> T(/+)+ X , (3.1) 

whereas Neutral Current (NC) interaction causes a modulation in the spectrum of the 
interacting neutrinos, 

vi{i?i) + N ^ vtfa) + X' . (3.2) 
These total cross sections can be written in terms of differential ones as follows 

= I Ev ^{E u ,y)dy , (3.3) 

o$,{K) = f Q ^{E v ,y)dy , (3.4) 

where E v is the energy of the incoming neutrino, mi is the mass of the outgoing charged 
lepton and y is the inelasticity parameter, defined as 

vccnc = i - "FT j ( 3 - 5 ) 

with Ei the energy of the outgoing charged (for CC) or neutral (for NC) lepton. 

lr The only exception is the resonant v e — > W~ production, occurring at _E Pc = 6.3 PeV, whose 
contribution to the total event rate remains nevertheless negligible |45|. 
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3.1 Deep inelastic neutrino cross sections 



-L/llcl j I vjrt: V J 


< Dec > 


< Unc > 


10 7 


0.2388 


0.2449 


10 8 


0.2180 


0.2223 


10 9 


0.2019 


0.2052 


10 10 


0.1900 


0.1928 


10 11 


0.1785 


0.1821 


10 12 


0.1542 


0.1601 



Table 1: Average inelasticity parameter for CC (< ycc >) and NC (< unc >) interaction 
for different incoming neutrino energy. 

Tau neutrino-nucleon cross sections have been calculated following the approach of 
Ref. [4*5] . based on the renormalization-group-improved parton model, and by using the 
most recent data on the parton structure functions of nucleons. 

The cross sections can be written in terms of the Bjorken scaling variables y and 
x = Q 2 /2MyE u , where — Q 2 is the invariant momentum transferred between the incoming 
neutrino and the outgoing lepton. Details of nucleon structure become important at very 
high-energy where available data are very poor or totally missing. As a consequence of 
this, the lack of knowledge of the parton structure functions at very low x (x <C 10~ 5 ) 
dominates the whole uncertainty on the cross section calculations at very high-energy. In 
the present analysis we have used the CTEQ6 [33] parton distribution functions in the 
DIS factorization scheme. 

The Q 2 -evolution is realized by the next-to-leading order Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi equations [H)]-[in]. The CTEQ6 distributions are particularly suitable for 
high energy calculations since the numerical evolution is provided for Q 2 in the range 10 0-8 
GeV 2 and for x down to 10 -6 (E v ~ 10 7 GeV). Values outside this range are calculated 
by extrapolation. 

The total cross sections for v T CC and NC inelastic scattering off an isoscalar nucleon, 
N = (n+p)/2 (n = # of neutrons, p = # of protons), are shown in Figure El The 
calculation for v T is not shown because antineutrino- and neutrino-nucleon cross sections 
become indistinguishable for E v > 10 6 GeV. This is due to the dominance at very high- 
energy, i.e small x, of the sea quarks on the valence ones. Figure 0] shows the distribution 
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Figure 3: The v T total cross section for CC (thick line) and NC (thin line) inelastic 
scattering off an isoscalar nucleon N = (n + p)/2 are here reported. The calculation has 
been performed by using the CTEQ6 Parton Distribution Functions (PDF) [33] , according 
to the prescription given in Ref. [13] . 

of the inelasticity parameter ycc = ^- — E T /E U , for different neutrino energies. An average 
inelasticity < y > can be defined by integrating the distributions given in Figure 0J Few 
relevant values of < y > are summarized in Table Q for CC (< ycc >) an d for NC 
(< Vnc >) interaction. 

Figure El shows a comparison between a CTEQ4-based parametrization of CC cross 
section and the corresponding calculation performed with CTEQ6. A substantial agree- 
ment is found up to 10 9 GeV, whereas a discrepancy of at most 30% at 10 12 GeV is 
observed (CTEQ4 prediction being larger). In this range of energy the uncertainty due 
to the lack of knowledge of parton distribution functions is expected to be overwhelming. 

4 The r energy losses 

Tau leptons with energy higher than 10 16 eV travelling in matter may loose a consistent 
fraction of their energy before decaying [20]. A precise knowledge of r energy loss is 
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Figure 4: Distribution of the inelasticity parameter, ycc = 1 — E T /E U , for different 
neutrino energies. 

therefore required in order to draw reliable predictions of expected signal rate at detec- 
tors. High-energy r's propagating through matter mainly interact by quasi-continuous 
(ionization) and discrete energy loss mechanisms, mainly by photonuclear interaction, 
direct electron-positron pair production and Bremsstrahlung. The ionization energy loss 
dominates for r-energy smaller than a few TeV, while radiative processes become relevant 
at higher energies. 

The direct electron pair production differential cross section has been calculated by 
Kelner and Kotov in the framework of QED theory [51]. We have used the well-known 
parametrization performed by Kokoulin and Petrukhin [52], which considers the cor- 
rections for atomic and nuclear form factors. The expression for the Bremsstrahlung 
differential cross section has been derived by Andreev and Bugaev [5H]- It takes into ac- 
count the nuclear structure of the target (elastic and inelastic form factors) and the exact 
contributions due to atomic electrons (screening effect and Bremsstrahlung on electrons). 

The complete formulas used here for computing electron pair production and Bremss- 



E >= 10 13 eV 
E v -10 14 eV 
E v -10 1B eV 
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Figure 5: The total i/ T -nucleon CC cross sections based on both CTEQ4 [UJ and 
CTEQ6 [33] are reported. 

trahlung cross sections are the ones reported in Ref. [51], with the only substitution of 
the muon mass with the r mass. Almost identical formulas are given in Ref. [55] and in 
Ref. [50] . where a slightly simplified formula for Bremsstrahlung (in agreement within a 
few percent with the one of Ref. [HI]) is actually used. A complete list of r matter cross 
sections written within the same theoretical framework adopted here is also given in Ref. 

The photonuclear differential cross section is calculated following the theoretical ap- 
proach developed in Ref. [S3]. According to this formalism, the cross section for photonu- 
clear interaction consists of two terms, where the first one, obtained within the Vector 
Meson Dominance Model, describes the non-perturbative contribution to the electro- 
magnetic structure functions. The parametrization given in Ref. [33] differs from the 
corresponding well known result of Ref. [HZ] by few new terms, negligible for muons but 
important for r's. The second term, with parameters updated with the most recent exper- 
imental data, describes the perturbative QCD contribution, not negligible at extremely 
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high energies (E T > 10 15 eV). The parametrization of these terms are provided up to 
10 9 GeV in Ref. jHEh whereas the values at higher energy are obtained by extrapolation. 
In order to derive the non-perturbative term we considered the recent accelerator data 
coming from the experiments ZEUS and HI [o^ IHO] . 

The average energy loss for a given discrete process k can be expressed in terms of the 
differential cross section, da k /dv, as follows 



where Na is the Avogadro's number, A is the mass number, v is the fraction of initial 
energy E T lost by the r at the occurrence of the process k and x is the thickness of the 
crossed matter, expressed in g/cm 2 . 

Figure El shows the (3 values for photonuclear interaction (the most relevant process 
at high energies), electron pair production and Bremsstrahlung versus r energy. Figure 
shows the energy loss due to individual electromagnetic processes and the total average 
energy loss defined as 



as a function of the initial, E T , and the final, E™ tn , r-energy. All results shown in 
Figures El -El are calculated for standard rock (Z=ll, A=22, g s = 2.65 g/cm 3 ). A detailed 
evaluation of the average (effective) range would require a full stochastic treatment of 
discrete interactions in order to correctly handle fluctuations. As shown in Ref. [SI] and 
[HQ] , fluctuations may strongly decrease the average range compared to the range of the 
average energy loss. Results of a full simulation of r propagation through matter are also 
shown in Ref. [HE! an d [62J . 

5 The Earth— skimming events 

Following the formalism developed in Ref. jM] , let be an isotropic flux of v T + V T . 
The differential flux of charged leptons emerging from the Earth surface with energy E T 




(4.6) 




(4.7) 



Finally, Figure |E1 shows the range of the average energy loss, 
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Figure 6: j3 value for photonuclear interaction (the most relevant process for r's at high 
energies), electron pair production and Bremsstrahlung. Results are shown for standard 
rock (Z=ll, A=22, g s = 2.65 g/cm 3 ). 



is given by 



d$ T {E T ,9,<j>) 



(5.9) 



dE T dtt J dE u dn 
where K(E V , 9; E T ) is the probability that an incoming neutrino crossing the Earth with 
energy E v and nadir angle 9 produces a lepton emerging with energy E T (see Figure Ell- 
in Eq. (j5.9j) . due to the very high energy of v T , we can assume that in the process 
v T + N t + X the charged lepton is produced along the neutrino direction. This 
process can occur if and only if the following conditions are fulfilled: 

a) the v T with energy E v has to survive along a distance z through the Earth; 

b) the neutrino converts into a r in the interval z, z + dz; 

c) the created lepton emerges from the Earth before decaying. 



a) The probability P a that a neutrino with energy E u crossing the Earth survives up a 
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Figure 7: The total average energy loss in standard rock is plotted (thick line). The 
contributions of each individual electromagnetic process are also reported (thin lines). 
Results are shown for standard rock (Z=ll, A=22, g s = 2.65 g/cm 3 ). 



certain distance z is 



where 



exp {-r ^ ' (5 ' 10) 



is the CC interaction length in the rock, g[r(9, z)\ is the Earth's density at distance r. 
The distance r is given by r 2 (9, z) = R 2 ^ + z 2 — 2R e zcos9, where R® ~ 6370km is the 
average Earth radius. The expression of P a does not take into account the atmosphere 
crossed by the v T before entering the Earth surface, because the CC interaction length in 
the air is almost three orders of magnitude larger than in the rock. 



b) The probability for v T — > r conversion in the interval [z, z + dz] is 

_ dz 



\ v cc (E u , 6; z) 



(5.12) 
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Figure 8: Range of average energy loss as a function of energy E T . E™ m has been fixed 
at 10 GeV. Results are shown for standard rock (Z=ll, A=22, g s = 2.65 g/cm 3 ). 



Here a comment is in turn. In order to produce a r emerging with enough energy to 
generate a electromagnetic shower detectable by FD, the charged lepton propagation in 
the rock is limited. Moreover, since an EeV neutrino has a \ v cc ~ 500 km, only quite 
horizontal v T will be able to produce detectable events. The charged current interaction 
will then take place near to the Earth surface where the average density is almost constant 
and equal to g s ~ 2.65 g/cm 3 . 



c) The probability P c that a charged lepton loosing energy survives as it travels through 
the Earth is described by the coupled differential equations: 



dPc 
dz 
dE T 

dz 



CT T E T 



Pc , 



= ~(Pt + lr E T ) E T S . 



(5.13) 
(5.14) 



Here m T = 1.8-10 9 eV, r T ~ 3.4- 10~ 13 s denotes the r mean lifetime, whereas the parame- 
ters P T ~ 0.71-10~ 6 cm 2 g _1 and 7 T ~ 0.35-10 -18 cm 2 g _1 GeV -1 , as discussed in Section 
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Figure 9: A neutrino u T crosses the Earth with energy E v at a nadir angle 9 and azimuth 
angle (p. Then it travels a distance z before converting into a charged lepton r, which 
emerges from the Earth surface with energy E T |31j . 



HI fairly describe the r energy loss in matter. The set of equations (j5.13|) . (j5.14|) can be 
solved by observing that, following the results presented in section El and shown in Table 
n the tau lepton produced at z carries an average energy which is a function of E v . Let 
£° = E°(E V ) = (1- < ycc >)E V be the transferred energy. By solving Eq.s ()5.13|) and 
(j5.14j) at the emerging point on the Earth surface one has: 

P c = (F(E v ,E T )f exp{ ^f^-^rl} , (5-15) 

{ CT T fj T Q s \E r E^(E U ) J J 

P T E*{E V ) exp{-g 8 /? T (2fleCos0-z)} 

T " /3 r + 7^ T (^)(l-exp{- & /3 T (2i? e cos^-z)}) ' l " J 

where 

F( ^'^ = E T (Pr + 7rE°(E„)) ' " = ^Pfg ■ (5 - 17) 

The above results improve the ones obtained in Ref. j^I], where a simpler parametrization 
for the t energy loss was adopted. 



The energy E T of the exiting lepton must be consistent with Eq. (|5.1fij) . This condition 
is enforced by the 5-function: 

p = J F _ [3rE° T (E u ) exp{-Q s (3 T (2R® C osd-z)} \ 
d \ T P T + It E^) (I -ex V {-g s (3 T (2Re cos 6 - z)}) ) ■ K ' ' 

By using the expressions for the different probabilities (|5.1U|) . (|5.12|) . (|5.15|l . and (|5.18l) 
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the kernel reads 

/■2i?0 cos 9 

K(E U , 6- E T ) = / PaPbPcPd dz . (5.19) 

Jo 

Once the integration over z is performed, one gets the simple result 

xexpj U^(^-—^]-2R @ cosea^(E v ) 6a N A ) , (5.20) 



ct t (3 tQs \E t E°(E v ) 



where 



/ a v r N r (E u )N A \ 

e = cc y - 1 • (5.2i) 

Eq. (|5.18J) requires that two conditions are fulfilled. The first one is E®(E U ) > J5 T (obvi- 
ously verified), while the second is 

cos9>cos9 min = — — \og{F{E v ,E T )) . (5.22) 

The expression (|5.20|) for the kernel K(E V , 9; E T ) leads to the total rate of upgoing r's 
showering on the Auger detector, and thus potentially detectable by the FD: 

at JE™ in JE* h Jco80 min aE v a\l 

x K(E„, 9- E T ) (l - exp {-^^}) e cos9d(cos9) , (5.23) 

where we have used the isotropy of the considered neutrino flux. In Eq. ()5.23j) the quantity 
S = 3000 km 2 is the geometrical area covered by the Auger apparatus, D ~ 10% is the 
duty cycle for fluorescence detection, E 1 ^ ~ 10 18 eV is the energy threshold for the 
fluorescence process, and E™ m is the minimum neutrino energy capable of producing a r 
at detection threshold. The quantity E™ ax is the endpoint of the neutrino flux. 

The exponential term in the r.h.s. of Eq. 1)5.23)1 accounts for the decay probability 
of a r (showering probability) in a distance H from the emerging point on the Earth 
surface. A detailed calculation of the fraction of Earth-skimming r's decaying inside 
the fiducial volume would require a full Monte Carlo simulation, which is not the aim 
of the present analysis. A reasonable way to estimate such a fraction is to require that 
r leptons travel a distance less than H before decaying. We choose for simplicity H = 
30 km, which is the radius of an emisphere approximatively containing the entire PAO 
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Figure 10: The effective aperture, A(E V ), as defined in Eq. (|5.25j ) is here plotted for Ef 1 = 
10 17 eV and H = 30 km, and for £f = 10 18 eV with H = 20, 30, and 40 km, respectively. 



apparatus. However, this ansatz provides a conservative estimate of the number of events. 
FD detection efficiency e has been taken by Ref. [S|. 

In Eq. (J5.23|) the integration over cos 6 can be easily performed and this yields to 



E rnax d& ( E ) 

D / dE u J* /* A{E V 



dt JE™ in dE u dVt 

where the effective aperture of the apparatus is here defined as 

7lS 



(5.24) 



A(E V ) 
x exp 



m T ( 1 



E t ((3 t + 1t E t ) 
I \ I f _ ( Hm T 



x 



uN 
CC 



ct t (3 t q s \E t E° t {E v ))\ 
\ + ^^\o g (F(E„,E T )) 



1 — exp 



ctv E T 



1 + 2 af Q S N A ) exp {-2 R e og G g s N A } {F{E„, E T )) a - 



(5.25) 



In Figure EH the quantity A(E U ) is plotted versus the neutrino energy for Ef 1 = 10 17 eV 
and H = 30 km, and for Ef = 10 18 eV with H = 20, 30, and 40 km, respectively. The 
aperture for E^ = 10 18 eV shows a maximum near E v ~ 10 19 eV, sensibly dependent on 
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the parameter H. For a lower energy threshold, for example 10 17 eV, A(E V ) would scale 
as shown in Figure ITU1 

A similar analysis for the PAO Surface Detector has been performed via a Monte Carlo 
simulation in Ref . [32] • in particular, one could compare the aperture given in Figure El 
with Figure 9 of Ref. [H2] ■ A direct comparison of the two calculations is difficult to be 
performed since the detector efficiency and the energy threshold of FD and SD differ. In 
our analysis we use an energy threshold E^ h = 10 18 eV. However, A(E V ) approaches the 
aperture of Ref. [32] as E^ decreases. 

In Table El the number of r-shower Earth-skimming events expected per year at the 
FD detector is reported for the models of neutrino flux discussed in section El 



dN T /dt at FD 


GZK-WB 


GZK-L 


GZK-H 


TD 


NH 


# of UP events/year 


0.02 


0.04 


0.09 


0.11 


0.25 



Table 2: Yearly rate of upgoing (UP) events at FD for different neutrino models given in 
Figures [T] and El 

In Ref. (HI] an explicit plot of the aperture as defined in our Eq. (|5.25J) is not given. 
However, it is worth observing that the number of events reported in Table EJfor TD model 
is in fair agreement with the corresponding result of Ref. once properly normalized to 
PAO effective area. Other cases are not immediately comparable since fluxes are different. 

In Figure ^2 the energy spectra of Earth-skimming events at FD are shown for the 
neutrino models of Figures H and El As expected, the maximum number of events is 
reached near to the assumed FD threshold (~ 1 EeV), even though the maxima of neutrino 
fluxes are at higher energy. This can be easily understood by observing that tau leptons 
emerging inside the Auger surface with large energy and almost horizontally will probably 
decay far from the apparatus and are thus undetectable by the FD. By using the integrand 
in the r.h.s. of Eq. (j5.24j) . the angle with respect to the horizontal of the Earth-skimming 
event a max = 9 max — n/2 can be determined. Here, 9 max denotes the nadir angle for which 
the kernel of Eq. (|5.24j) has the maximum. This quantity is a function of neutrino energy 
via Oqq and is shown in Figure fT2l 

As pointed out for example in Ref. [2E] , v T crossing deeply the Earth could experience a 

17 




regeneration phenomenon which would eventually let r's propagate up to the detector(see 
Figure IT^j). Furthermore, this process would also lead to the production at the CC vertex 
of secondary electron or muon neutrinos, then decaying in electrons or muons. 

The contribution of regenerated z/ T 's to the Earth-skimming event rate is negligible, 
since the process shown in Figure Ell is a second order transition in the weak coupling 
constant. Moreover, for Ef 1 > 5-10 17 eV the Earth strongly suppresses the expected signal 
due to t energy losses [2SH2HI- As shown in Ref. [23123, the value of E 1 ^ plays a crucial 
role in the estimate of the regenerated neutrino event rate, which results negligible for 
threshold energies larger than 10 18 eV. Higher order processes including NC interactions, 
would also give negligible modifications to the numbers presented in Table 121 Finally, 
the associated secondary electron and muon neutrino flux does not enhance the Earth- 
skimming signal at PAO, since the charged leptons they produce are either absorbed or 
not detectable by FD. 

As well known, neutrino induced extensive air showers (see Ref. jBB] for an extended 
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Figure 12: The most likely exit angle with respect to the horizontal for the emerging r is 
plotted versus the energy of the primary neutrino. 




Figure 13: Second order regeneration process for v T . 
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discussion) can be disentangled by the ordinary cosmic ray background only for very 
inclined and deep showers. An estimate of the expected number of downgoing events 
(DW) within 30 km from the FD detector is given in Table E3 for different fluxes and 
increasing zenith angle. The DW event rates result to be comparable with the Earth- 
skimming ones of Table El As a final remark, the DW and Earth-skimming event rates 
are both affected by the uncertainty on the z/-nucleon cross section, but, unlike the DW 
events, the Earth- skimming rate is also affected by the uncertainty on the tau energy 
losses in rock. 





GZK-WB 


GZK-L 


GZK-H 


TD 


NH 


# of DW events/year (> 60°) 


0.04 


0.08 


0.22 


0.25 


0.54 


# of DW events/year (> 70°) 


0.02 


0.04 


0.10 


0.12 


0.25 


# of DW events/year (> 80°) 


0.005 


0.01 


0.03 


0.03 


0.07 



Table 3: In the table are reported the expected number of downgoing events (DW) within 
30 km from the FD detector for different fluxes and increasing zenith angle. 

6 Conclusions 

Ultra High Energy i/ r 's [E v > 1 EeV) could have real detection chance at giant sur- 
face apparatus like PAO. Almost horizontal tau neutrinos crossing distance in the Earth 
of the order of their interaction length might produce r-shower Earth-skimming events 
potentially detectable by the Auger FD. 

In this paper for a representative sample of u T fluxes either produced as cosmogenic 
neutrinos or in more exotic scenarios, the number of r-shower Earth-skimming events 
per year expected at the FD of the Pierre Auger Observatory has been computed. For 
this calculation we have used the CTEQ6 [33] parton distribution functions in the DIS 
factorization scheme and recent estimate of radiative r energy losses in matter. A decrease 
of 

a cc allows r to emerge with a smaller nadir angle (less horizontal) and thus increases 
the number of Earth-skimming events. In the relevant region of E u , we essentially have 
dN T /dt oc (o^cO -1 ; thus a factor 1/2 in the neutrino- nucleon cross section leads to an 
almost double number of events at FD. As a further remark, we point out that the level of 
theoretical uncertainties on cross-sections at extremely high energy is very large (like in 
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case of new physics above the TeV energy region), due to the poor experimental knowledge 
of parton density functions for very small x. 

The number of events per year at the FD of the Pierre Auger Observatory is presented 
in Table El for the different neutrino models reported in Figure ^ and El For a five years 
PAO detection campaign (South plus North) we essentially have to scale these numbers by 
a factor ten and thus at least the exotic models could produce detectable events. However, 
it is worth noticing that this estimate has been performed under a rather conservative 
assumption on the possibility of r showering on the apparatus: we have used an average 
value for the parameter H which essentially allows only one particle out of three to 
decay inside the detection fiducial volume. In a paper in progress a more sophisticated 
simulation is going to be performed by taking into account the real morphology of the 
PAO site. This analysis should slightly increase the numbers presented in Table which 
however give a reasonable estimate of the number of events at FD. 
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